SCIPP-12/04 



Diphotons from Tetraphotons in the Decay of a 125 GeV Higgs at the LHC 

Patrick Draper 1 ^ and David McKeen 2 '[j] 

1 Santa Cruz Inst, for Particle Physics, Univ. of California, Santa Cruz, CA 95064, USA 
2 Department of Physics and Astronomy, University of Victoria, Victoria, BC V8P 5C2, Canada 

Recently the ATLAS and CMS experiments have presented data hinting at the presence of a 
Higgs boson at mj, ~ 125 GeV. The best-fit h — > 77 rate averaged over the two experiments is 
approximately 2.1 ± 0.5 times the Standard Model prediction. We study the possibility that the 
excess relative to the Standard Model is due to h — > aa decays, where a is a light pseudoscalar that 
decays predominantly into 77. Although this process yields 47 final states, if the pseudoscalar has 
a mass of the order tens of MeV, the two photons from each a decay can be so highly collimated 
that they may be identified as a single photon. Some fraction of the events then contribute to an 
effective h — > 77 signal. We study the constraints on the parameter space where the net h — > 77 
rate is enhanced over the Standard Model by this mechanism and describe some simple models that 
give rise to the pseudoscalar-photon interaction. Further tests and prospects for searches in the near 
future are discussed. 
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I. INTRODUCTION 

Determining the nature of electroweak symmetry 
breaking (EWSB) has been a primary goal of parti- 
cle physics for several decades. In the Standard Model 
(SM), EWSB occurs when the neutral component of a 
single scalar weak isospin doublet possesses a vacuum 
expectation value (VEV). The Higgs boson corresponds 
to the physical excitations of this field. Recent ex- 
perimental advances make it clear that we are enter- 
ing a new phase in searches for the Higgs boson, with 
tantalizing excesses appearing in several SM-like Higgs 
search channels at both the CERN Large Hadron Col- 
lider (LHC) and the Fermilab Tevatron. These hints 
point to a relatively light state with a mass in the range 
122.5 < m h < 127.5 GeV [TJ|2]. In the SM, such a light 
Higgs boson has an extremely narrow width of the order 
10~ 5 x nih- The tiny width of the SM Higgs makes it a 
sensitive probe of new physics (NP) beyond the SM, es- 
pecially in the only loosely constrained scalar sector [3J. 
New states coupled to a SM-like Higgs can have appre- 
ciable effect on its decays without ruining the excellent 
agreement of the SM with data observed elsewhere thus 
far (see, e.g., [1]). 

In fact, the excesses seen at the LHC and Tevatron may 
already be pointing toward interesting deviations from 
the SM predictions for the Higgs branching ratios. The 
largest statistical power thus far comes from the searches 
for h 77 at the LHC, where the ATLAS experiment 
finds a 2.8a excess at rr%h ~ 126 GeV with a best-fit 
signal strength relative to the SM of approximately 2.0 ± 
0.8 [U [5]. The CMS experiment finds a 3.1cr excess in 
h —> 77 at rrih ~ 124 GeV with a best-fit signal strength 
relative to the SM of approximately 2.1±0.6 Under 
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the assumption that these excesses are due to the same 
new particle, we can naively combine the diphoton rates 
from the two experiments and estimate 

Rate(/i^7 7 )/SM~2.1±0.5. (1) 

This estimate is not rigorous and the error bar is far 
from conclusive. However, the uncertainties will continue 
to decrease as data accumulate, and the central value 
may very well remain high. Therefore, it is of interest 
to classify models that can alter the h — > 77 signal and 
study their other predictions. A number of recent papers 
have contributed to this program, including discussions of 
the effects of superpartners on the diphoton rate [7] , the 
effects of more general new fermion and scalar states [8] , 
the effects of singlet-doublet mixing 9], the predictions 
in the case of minimal universal extra dimensions |10) . 
interference effects from charged Higgs contributions [TT] , 
and the possibility that the signal is due to the Randall- 
Sundrum radion |12j . to name a few. 

In this paper we consider a different scenario in which 
a 125 GeV Higgs boson can appear to have a larger 
branching ratio into 77. We introduce a very light pseu- 
doscalar of mass in the range m a E (~ 10 MeV, ~ rn„), 
and we allow the Higgs boson to decay in the channel 
h — > aa. The pseudoscalars produced in these decays are 
extremely boosted in the lab frame, causing their decay 
products to be quite collimated. Furthermore, for such 
light pseudoscalars, there are no kinematically available 
hadronic decay modes. As first studied in [13], the de- 
cay a — > 77 induced by the effective coupling aF^F^ 
can easily dominate, leading to a h — > aa — > 47 signa- 
ture (in contrast to the "buried Higgs" scenario where 
the dominant decays are hadronic [14] ). Because of the 
large boost for the pseudoscalars produced in h decays, 
the photon pairs from each a can be so highly collimated 
that a significant fraction may be identified as single pho- 
tons, even in the ATLAS and CMS detectors which are 
very good at distinguishing two closely separated photons 
(e.g. originating from ir° or r\ decays). 
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A related scenario, although not in the context of Higgs 
boson decays, was considered in [15], which studied de- 
cays of a new heavy vector through pseudoscalars to ex- 
tremely collimated "photon jets." Since their vector mass 
is much larger than 125 GeV, they consider much heavier 
pseudoscalars. Our analysis of the conditions necessary 
for photon jets to be identified as single photons is com- 
plementary to that of |15) and provides further motiva- 
tion for the detailed study of photon jets. Previously, the 
h — > aa —> 47 scenario was considered in [4j [13l [16] and, 
in particular, the relevance of this channel to h — > 77 
searches at the Tevatron was studied in [T3]. While 
Ref. [13] focused mainly on pseudoscalars heavier than 
200 MeV and heavier Higgs bosons, we will find that a 
significant Aj — > 2j fake rate at the LHC with a 125 GeV 
Higgs boson requires a pseudoscalar lighter than about 
the mass of the ir°. In contrast to previous works, we also 
study the interplay between the net expected h — ¥ 77 sig- 
nal and the signal in other Higgs search channels, provide 
an analysis of the compatibility of the light pseudoscalar 
scenario with the current Higgs data from the LHC and 
Tevatron, and study in detail the constraints from LEP 
and low-energy experiments. For further discussion of 
Higgs decays to light pseudoscalars, see, e.g., [T7] . 

This paper is organized as follows. In SecjHjwe discuss 
the basic features of the model. In Sec. |III| we analyze 
the changes in the expected rates for SM Higgs search 
channels at the LHC. We give a detailed estimate of the 
probability for a collimated photon pair to be identified as 
a single photon, and perform a basic statistical analysis of 
the constraints on the light pseudoscalar parameter space 
coming from current Higgs searches at the Tevatron and 
LHC. In Sec. |IV| we study the constraints following from 
the muon anomalous magnetic moment, direct searches 
at LEP for e + e~ — > 7 + inv., meson decays, beam dump 
experiments, and nuclear physics. In Sec. [V] we survey 
model building for the light a and its coupling to photons, 
and consider the possibility that the latt er is generated 
by r loops, or by a — ir° mixing. In Sec. VI we discuss 
the prospects for directly producing pseudoscalars in the 
mass and coupling range of interest at Primakoff-type 
experiments. In Sec. |VII] we conclude, and an Appendix 
generalizes our main results to the case where a has a 
substantial invisible branching fraction. 



II. MODEL 



has a small mass m a coming from small explicit breaking 
of the symmetry, and it interacts with the SM via 
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/: iul = —^afrfH-^aF^F^, (2) 



where H is the SM Higgs doublet, is the photon 
field strength, and e is the positron charge. A and M are 
scales describing the strength of the higher dimensional 
operators. 

Expanding around the Higgs VEV v ~ 246 GeV, H° = 
(v + h) I a/2, leads to a rate for the Higgs to decay to 
pseudoscalars of 
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where we have assumed that m a <C m^. Using Eq. |2|) 
we can also calculate the rate for a to decay to 77, 



r (a -> 77) 
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= 2.68 x 10~ 8 MeV 
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Pseudoscalars produced by the decay of 125 GeV Higgs 
at rest then have a decay length of 



7CT ~ 1.15 cm 



M 



10 GeV 



)'( 



m„ 



40 MeV 
m h 



(5) 
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If we require > 90% of the decays to occur inside the 
electromagnetic calorimeters of ATLAS and CMS (~ lm 
radius) , jct should be less than about a half meter. We 
can express the scale M in terms of the decay length, m a , 
and m/j, 



M = 9.3 GeV 



(S" 2 ( 



40 MeV 

m h 



(6) 
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125 GeV 
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In the Appendix, we expand this simple model to al- 
low for a to additionally decay invisibly. We discuss the 
benefits and further constraints that this entails there. 



In addition to the SM, we consider a real pseudoscalar 
a which is the pseudo Nambu-Goldstone boson of some 
spontaneously broken approximate global symmetry^] It 



1 A frequently-discussed model containing a light PNGB is the 
NMSSM in the approximate PQ- or R-symmetric limits. In that 
model the pseudoscalar has sufficiently large fermionic couplings 
that in the mass range we consider, it is completely ruled out by 
low-energy experiments such as those considered in Sec. II VI f 8 . 



III. HIGGS SIGNAL 

The interactions described in Sec. [H] do not apprecia- 
bly affect the production cross section of the SM Higgs 
boson at the LHC. In this section we study the effects of 
the decays h — > aa, a — > 77 on the expected Higgs-to- 
diphotons rate. We also estimate the constraints on such 
decays from existing SM Higgs searches at the Tevatron 
and LHC. 
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A. Modified Rates in SM Higgs Search Channels 

We define the ratios of the Higgs branchings to photons 
and to SM fermions, /, or gauge bosons, V = W, Z, to 
their values in the SM as 

B (h -> 7 7 ) cff = i? 77 x B S m (h -> 77) , (7) 
B{h-t ff, VV) = R xx x B SM (h -> //, W) . (8) 

The presence of the h — > aa decay suppresses uniformly 
the rates into all SM final states, yielding 

Rxx = l-B(h-^aa). (9) 

Note that to avoid cluttering notation, and because the 
rescaling is the same in both cases, we use Rxx to rep- 
resent both the fermionic and W, Z rescalings relative to 
the SM. 

The suppression in Eq. [9] is also present for the pure 
77 final state; however, the subscript "eff" in Eq. Q 
indicates that in pseudoscalar models there can be an 
additional, effective contribution to the measured h — > 77 
rate. Since the LHC detectors have finite resolution, a 
fraction of h — » aa — » 47 events may have sufficiently 
boosted photon pairs that each pair is identified as a 
single photon. Assuming a SM production cross section 
for h and a 100% branching of a — > 77, the measured 
diphoton branching ratio will be 

B (h -> 7 7 ) cff = B (h -> 77) + e x B (h -> aa) (10) 

or 

^ 7 = l + g (/w aa) ( gsM( ^ 77) -l). (11) 

In these formulae, e is the probability that both 77 pairs 
are identified as single photons, so that four photons ap- 
pear as two. We see that to achieve an effective diphoton 
rate greater than or equal to the SM rate requires 

e > B S m (h -> 77) ^ 0.0023 (12) 

for m h = 125 GeV. 

B. 47 — > 27 Misidentification Rate 

Estimating e is complicated by several factors. We 
would like to be conservative in our estimate of the ex- 
pected rate; on the other hand, underestimating the rate 
by too much might falsely indicate that some model 
points are allowed, when in fact the true rates are so 
large that the points are already ruled out by the LHC. 
We base our estimate of e on the ATLAS selection cri- 
teria used to identify isolated photons in their cut-based 
analysis, and we comment on differences with CMS. 

ATLAS uses a number of calorimeter variables to 
parametrize the shape of an electromagnetic shower, 



which can then be used to discriminate true isolated pho- 
tons from backgrounds. The background from isolated 
7]-° — >. decays bears strong similarity to our a — > 77 
process, and ATLAS efficiently vetoes isolated pions us- 
ing information from the first layer of the calorimeter, 
which has finely-segmented strips in 77 |19j . The primary 
discrimination variables in the first layer are E ratio , which 
is the difference in energies between two strips containing 
energy maxima normalized to their sum; AE, which mea- 
sures the difference in energies between the strip with the 
second-largest energy maximum and the strip with the 
minimum energy between the first two maxima; F S idc, 
which is the fraction of the energy deposited in seven 
strips in 77 around the maximum that does not fall into 
the central three strips; w s3 , which measures the energy 
deposition in the two strips adjacent in 77 to a strip with 
an energy peak, relative to the total energy in the three 
strips; and w sto t, which generalizes w S 3 to approximately 
twenty strips in 77 and two strips in <fr |20j . 

To simplify our analysis, we begin by restricting our 
attention to photons that do not convert to e + e~ pairs 
in the tracker, and which are so highly collimated that 
a second energy maximum does not appear in the first- 
layer calorimeter strips. In this case AE is set to 0, -E ra tio 
is set to 1, and most of the energy will be deposited into 
just a few adjacent strips. Therefore we expect that the 
most sensitive variable will be w S 3, defined precisely as 



w s3 = j^E i {i-i max Y/Y J E i , (13) 

V i i 

where i labels the strips. On average the photon pair 
from a ir° decay generates a larger w S 3 value than a true 
single-photon event and thus may be efficiently rejected. 
On the other hand, the probability that a photon pair 
passes the w S 3 cut should increase with the collimation 
of the pair. To easily estimate e, we would like to ap- 
proximate the cut on w S 3 by a cut on the photon pair 
separation. 

For unconverted events, ATLAS uses a weakly 77- 
dependent cut on w s3 that is typically between 0.6 — 0.7, 
and is about 0.66 for the most central strips in the bar- 
rel. The average value of w S 3 for true photons in these 
strips is about 0.58 [21]. We can reproduce this num- 
ber with the following simple model. First, we assume 
that a single photon lands in the center of one of these 
most central strips, and deposits its energy according to 
a Gaussian distribution in 77 with standard deviation 0.52 
times the smallest strip width. This value is chosen so 
that the photon deposits about 70% of the energy into 
the central strip and 15% into each adjacent strip, giving 
a Ws3 value of 0.58. Subsequently, we compute w S 3 for 
a pair of photons as a function of 77 separation, averag- 
ing over the impact point in the central strip. We find 
that for A?7 77 = 0.0015, 77j s3 ~ 0.66, equal to the ATLAS 
cut in these strips. On the other hand, the central strips 
have a width in 77 of A77 str ip = 0.0031 (the smallest in 
the calorimeter.) Therefore, we estimate that the cut on 
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u> S 3 can be simulated by a cut on the photon separation, 
given by 



A?7 77 < 1/2 x A?7 strip . 



(14) 



For larger 77, where larger strips are present, we also re- 
quire a separation less than 1/2 of the relevant strip size, 
which is probably mildly conservative since on average 
the energy leakage into adjacent strips from a single- 
photon hit will be less. 

There are two simple ways in which w S 3 may become 
insufficient. First, photons that are less highly collimated 
(such as those that appear in the decays of more massive 
pseudoscalars) may deposit energy in strips that are suffi- 
ciently separated (more than three strips apart) that w S 3 
becomes insensitive. We assume that such events are very 
efficiently rejected by the other first-layer discriminators 
described above, and therefore our cut on A?7 77 is still 
a good proxy for the cuts on those variables. Secondly, 
photons may be closely spaced in 77, but more broadly 
spaced in </>. Since the segmentation in <p is much coarser 
than in 77, such events must be very separated in <fi in 
order for cuts on the variable (which measures the 
energy distribution across several second-layer (f> cells) to 
reject them. A cut on A</> 77 can simulate the cut on R^, 
but since such events are geometrically rare, the net effi- 
ciency is quite insensitive to the precise value of this cut. 
For dcfinitencss we set the cut on A0 77 to be equal to 
the smallest second-layer cell size in <j> as a function of 77 
(A^ ce ii = 0.025 for 77 < 1.4): 



A0 77 < A<j> c 



(15) 



So far we have discussed only events which contain 
no converted photons. Conversions happen with an 77- 
and -&r-dependent probability that ranges from about 
10% at low 77 to more than 50% at larger 7/ [20.. Since 
h — > aa events produce four photons in the final state in- 
stead of two, a larger fraction of our events will contain 
at least one converted photon than are contained in pure 
h — > 77 events. In our conversion events, one cluster 
in the calorimeter may contain one e + e~ pair and one 
7, or two e + e~ pairs, so we might imagine that these 
events would be easy to distinguish from pure single-7 
events and may even be vetoed. In the case with 7e + e _ , 
there is a mismatch between the momentum measured by 
the tracker (sensitive only to the charged particles) and 
the energy deposit in the calorimeter. In the case with 
2e + e~, two conversion vertices may be reconstructed in 
the tracker. Both of these signatures have been con- 
sidered for rejecting isolated tt° backgrounds [55]; how- 
ever, at present, the ATLAS photon identification analy- 
sis does not use an E/p cut, and does not automatically 
veto events with multiple conversion vertices |20) (only 
the "best" reconstructed conversion vertex is used, and 
is determined by the conversion radius and the number 
of tracks associated with the vertex.) Furthermore, the 
cuts on the calorimeter variables are relaxed somewhat 
to accommodate the fact that the energy deposits for 
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FIG. 1. Fraction of events where both photon pairs are suffi- 
ciently collimated to pass ATLAS isolation cuts, as a function 
of the pseudoscalar mass. 



single-photon conversions tend to spread mildly in 77 and 
considerably in <f> (due to the magnetic field.) Therefore, 
we will make the approximation that the value of e rele- 
vant for 47 events containing conversions is the same as 
the value of e for the unconverted sample; namely, it is de- 
termined by applying the collimation cuts in Eqs. ( 14|15 1 
to the sample of parent photons before they convert. 

In summary, for all classes of photon events, we make 
the approximation 



e - e C oll, (16) 

where e co n is the rate at which the 47 events are expected 
to pass the 77-dependent collimation cuts At7 77 < 1/2 x 

A?7 s trip, A</> 77 < A0coii. 

The collimation of a photon pair produced in an a 
decay is controlled by the ratio 2m a /E a , where E a is 
determined by the momentum of the Higgs boson. We 
simulate 7 TeV gg — >■ h — >■ aa —> 47 events in Madgraph 
5 [53]) an d on the subsample of events with all photons 
satisfying \rj\ G (0, 1.37)U(1.52, 2.37) (the region in which 
ATLAS reconstructs photon candidates for the h — > 77 
analysis) , and which pass the ATLAS photon px cuts [5] , 
we compute the efficiency e co ii- In Fig. [T] we plot e co n as 
a function of 7n a j^] 



2 Production mechanisms other than gluon fusion result in some- 
what different kinematic distributions for h. However, gluon fu- 
sion dominates the production with O(10%) contribution from 
all other channels [5], so the effects on e co n from these processes 
will be second-order. We neglect them in this study. In addi- 
tion, we have checked that the effect of raising y/s to 8 TeV is 
negligible. 
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To cross-check that our cuts effectively reproduce the 
rejection power of w s ^ and the other discrimination vari- 
ables, we use the known isolated ir° rejection power of 
ATLAS as a function of (Et,t]) (given in Fig. 11 of 
Sec. 5 of |19| ) to estimate e from our simulation at the 
mass point m a = m,o. We find good agreement be- 
tween the e derived from this method and e co ii obtained 
from the collimation cuts. Since the ATLAS tt° rejection 
power was defined on events containing both converted 
and unconverted photons, we are reassured that Eq. ( 16 ) 
is plausible. 

Unlike the ATLAS analysis, the CMS analysis does 
place cuts on the ratio of the calorimeter energy mea- 
surement to the momentum measured in the tracker in 
order to isolate single photons from ir° decays [21]. The 
cut may reduce e relative to e co n- However, the CMS 
ECAL has a barrel granularity ~ 6x larger in r\ than 
that of ATLAS [25], increasing e co n. For this analysis we 
assume that the net efficiencies e will be comparable be- 
tween the two detectors, but in principle any difference 
between them translates into an experiment-dependent 
and possibly conversion-dependent prediction for the 77 
rate and could be used to probe the value of m a . 

To summarize this section: 

(1) High collimation is sufficient for a significant num- 
ber of 27 clusters to be misidentified as single isolated 
photons by the experimental analyses. 

(2) In order to produce an interesting misidentifica- 
tion rate, m a < m^o. For larger m a there may still be 
misidentified events in the h — > 77 signal, particularly if 
B(h — > aa) is not small, but the scenario should be more 
properly studied in the context of an h — > 47 signal. 

(3) Contamination of h — > 77 by misidentified 47 
events will be visible as an excess of the total number 
of conversion events relative to the 77 expectation, as an 
excess of conversion events with multiple reconstructed 
vertices, and as an excess of conversion events with a 
mismatch between the track pj< and the energy deposit 
in the ECAL. 



C. Current Constraints and Fits 

In this study we are interested in the possibility that 
the current excesses in Higgs searches at ATLAS, CMS, 
and the Tevatron are due to a Higgs boson with mass 
near 125 GeV. Under this hypothesis, some parts of the 
(rn a ,B(h aa)) parameter space are disfavored: they 
produce too few vector boson or fermion events, or too 
many or too few diphoton-like events. We estimate the 
constraints with a matched-filter, taking each point in 
parameter space as a template and estimating a best-fit 
amplitude R for the signal strength at the point In the 
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FIG. 2. Regions of light pseudoscalar parameter space that 
are favored (blue/light) and disfavored (red/dark) by the cur- 
rent best-fit signal strengths in the h — > 77, ZZ, WW, bb, tt 
channels. Contours are overlaid for the net diphoton (solid 
green) and ZZ, WW, bb, tt rates (dashed yellow) expected at 
the LHC relative to the SM rates. 



Gaussian limit R can be estimated as 
R = a\Cr/d,A, 



(17) 



where dj is the "data," which we take to be the best- 
fit amplitudes relative to the SM for the channels h — > 
77 [5], h ZZ 41 123, h ->■ WW -> Ivlv [28], 
{W/Z)h (11, Iv, vv)bb [25], and h^t tt [3D], presented 
by ATLAS at m h = 126 GeV [3T] , the channels gg -5> 
h -> 77 [32], qqh 9977 [i |32], h -> ZZ _-> 41 [33], 
h -> WW -> Ivlv [34], (W/Z)h -> (ll,lv,vv)bb [35], and 
h tt [3B], presented by CMS at m h = 124 GeV [37], 
and the channels h — > bb and h — > WW presented by 
CDF and DZero at m h = 125 GeV [38]. For ATLAS 
and CMS we chose slightly different mass points based 
on where their respective 77 excesses are most signifi- 
cant; we assume that the current experimental resolu- 
tion is large enough that both excesses can come from 
the same Higgs-like particle. The vector i$ gives the 
theoretical prediction for each data point (either i? 77 or 
Rxx), and is the inverse covariance matrix set by 
the squared symmetrized error bars (we include also cor- 
relations which we estimate from the luminosity uncer- 
tainty and the theoretical uncertainty on the gluon fusion 
cross section). The error on the estimate R is given by 



a^iUC-M,)-^. 



(18) 



3 For other recent studies of the Higgs best-fit cross section data 
using similar \ 2 analyses, see |26| . 



Model points where R — 1 is outside the 90% CL band 
around R are then disfavored. Of course, this statistical 
procedure is quite approximate at this stage, and serves 



only to get an estimated exclusion region if the data sets 
are large and nth actually is at 125 GeV. A precise cal- 
culation would construct the full Poisson likelihood with 
all correlations. With these caveats in mind, in Fig. [2] we 
present the regions of (m a ,B(h — > aa)) parameter space 
disfavored by the current LHC and Tevatron results, in- 
cluding contours of i? 77 and Rxx- The contours make 
clear the fact that the disfavored regions are mainly con- 
trolled by the excesses in the diphoton channels: increas- 
ing it above a few times the SM rate is in tension with the 
current excesses. We also maximize the likelihood over 
the plane, fixing the signal strength parameter to R = 1 
and constraining m a < m % . The best-fit point lies near 
the i? 77 = 2 contour at the m w boundary; relaxing the 
constraint, it would move beyond to near the intersection 
of the i? 77 = 2 and Rxx = 0.5 contours. However, the 
data are over-fit, x 2 /d.o.f. = 7.3/10, and the statistic 
is very shallow, particularly in the direction of constant 
i? 77 . Therefore, the precise location of the best- fit point 
has little meaning, and we do not show it in Fig. [2j For 
comparison, we find that the SM x 2 /d.o.f. = 12.2/12. 
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FIG. 3. A representative diagram of the leading contribution 
of the pseudoscalar, a, to (g — 2) M . 

stronger limits on the parameter space of the model and 
comment that these limits can likely be circumvented 
with relatively straightforward extensions of the model 
under consideration. 



A. Muon Anomalous Magnetic Moment 



IV. DIRECT CONSTRAINTS ON a 

We now discuss existing experimental constraints on 
the model in Eq. ([2]). Similar considerations have been 
undertaken with light pseudoscalars in [TS] [3<|] . 

Constraints from e + e~ — > 7a, quarkonia decays, and 
beam dump experiments are robust in the sense that 
they are controlled largely or entirely by the coupling 
of a to photons. Constraints from the muon anomalous 
magnetic moment and meson decays involving flavor- 
changing neutral currents are sensitive not only to the 
aFF coupling, but also any small coupling of a to SM 
fermions that might be present in the Lagrangian. These 
contributions may interfere constructively or destruc- 
tively with the aFF terms- we can even envision sit- 
uations where some constraints are mitigated by finely 
tuned tree level couplings of a to SM fermions. We as- 
sume for the purpose of this section that the fermionic 
couplings at the scale M are small enough that the lead- 
ing contributions to the constraints come from photon 
loops. 

Since aFF is dimcnsion-5, the amplitudes are typically 
logarithmically divergent. To remove the divergence, we 
cut the integrals off at momenta ~ M , which is moti- 
vated by the view that the aFF interaction is an effective 
coupling resulting from some new physics at a scale M — 
whether a is a composite with constituents with masses 
of that order or is fundamental and coupled to photons 
via particles with such masses. The phenomenological 
need for M to be tens of GeV poses some model-building 
puzzles in this respect and we will sketch a few potential 
scenarios leading to such a scale in Sec. |Vj 

Taken together, these low energy constraints should 
be viewed mainly illustratively. We attempt to be as 
conservative as possible, erring on the side of presenting 



A light pseudoscalar that couples to two photons will 
contribute to the anomalous magnetic moment of the 
muon. The leading correction, shown in Fig. [3j occurs 
at order a 3 and is analogous to the 7r°-pole portion of 
the hadronic light-by-light contribution to (g — 2) M . 

To estimate this effect as a function of m a and M, we 
make use of the expression for the 7r°-pole contribution 
to (g — 2)^ in [40] , rescale it by the strength of the photon 

coupling relative to that of 7r°, (47r 2 F 7r /Af) 2 , replace the 
pion mass with m ai and cut off the logarithmic divergence 
at the scale M instead of at ~ ATT 2 F n . To set limits 
we require that this contribution is less than the current 
deviation between experiment and theory on Aa M = (g — 
2)^/2 of about 25 x 1CT 10 gj. This sets a lower limit on 
M around 47r 2 F 7r = 3.64 GeV that is not very sensitive 
to m a for the range of pseudoscalar masses we consider. 

Because of the lepton mass dependence of corrections 
to g — 2, the electron anomalous magnetic moment is less 
constraining and does not appear in our limits. 

B. e + e~ — > Limits 

The reaction e + e~ — > ja can proceed through the s- 
channel exchange of a virtual photon. The cross section 
for this process is independent of the center-of-mass en- 
ergy, 

da _ 2ir 2 a 3 3(1 + cos 2 0) 
dcos6~ 3M 2 8 ' (1 ' 

where 9 is the angle between the photon and the beam 
axis in the center-of-mass. If a decays without being 
detected, then this process is subject to limits from 
e + er — > 7 + inv. We apply limits from the DELPHI 
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FIG. 4. Diagram that gives the leading contribution to s — > 
d + a from an effective interaction between a and the top 
quark. 

collaboration using 650 pb^ 1 of data at center-of- 
mass energies y/s = 180 GeV-209 GeV|^] In the photon 
energy range of interest for m a -C -^/s, E 1 ~ \/s/2, we 
require that there are fewer than 20 events and assume 
that the product of angular acceptance and efficiency 
is ~ 0.5 . Assuming that the a decay length must be 
larger than 2 m to be unseen at DELPHI and a boost of 
7 ~ 90 GeV/TO a excludes the following region of param- 
eter space: 

— < M < 125 GeV. (20) 

40 MeV/ ~ v ' 



C. Meson Decays 

The interactions in Eq. |2| lead to effective couplings 
of a to SM fermions at the one-loop level. This can have 
effects on rare meson decays involving photons or missing 
energy if a decays too late to be detected. 

Effective couplings to up-type quarks can then lead 
to flavor-changing transitions in the down sector such 
as s —> d + a as in Fig. [4] This leads to the flavor- 
changing decay K — > 7r ± a which must be confronted 
with experimental data on K — > tt ± jj. The amplitude 
for s — > d + a can be very roughly estimated as 

, x a 2 , 2 fM 2 \ [G F m 2 \ 
M{S ^ d + a) ~M l ° S (21) 
x V tdVts (m s d L s R - m d d R s L ) , 

where we have again cut off the logarithmic divergence at 
the scale M. Using this we arrive at a rate for K ± — > 7r ± a 



4 We note here that limit on e + e — — > 7+inv. from the ASP experi- 
ment 43 that has been used in previous studies to constrain light 
pseudoscalars coupled to photons do not actually apply. This is 
because their analysis imposed a cut on the photon's energy of 
E-y < y/s/2 — 4.5 GeV = 10 GeV to eliminate backgrounds from 
e+e — — > 77, rendering it insensitive to invisibly decaying parti- 
cles produced in association with a photon with a mass less than 
16 GeV. 



of 

r [K^ -> 7t eti) ^ bg 4 ^ £ * (22) 

V ' M 2 \m 2 J \ 2 10 tt 5 J V ; 

mr K \ m s - m d ) 1 1 

where the form factor /g (0) ~ I [Hj and is the pion 
momentum in the kaon rest frame. We apply the limits 
on hypothetical particles with masses less than 100 MeV 
decaying to two photons from [JS] (generally at the 10~ 7 
level) to B (K+ -> ir+a), resulting in M > 1-2 GeV. In 
addition, if a is sufficiently long-lived so that its decay 
is not detected, it will give a contribution to the process 
K + — > tt + + inv. The current experimental upper limit 
on the branching ratio for K + to decay to ir + and a light 
invisible particle is 7.3 x 10~ n [46] . For this limit to 
apply, a must have a decay length larger than about 1 m 
when produced with a 7r + in the decay of a K + at rest. 
The excluded region is then 

1500 GeV f — <M< 17 GeV. (23) 

V40 MeV/ ~ y J 

Analogous limits from B ± — > K ± transitions are less 
stringent in the pseudoscalar mass range of interest. 

Si quarkonium states can decay to 7a through an s- 
channel virtual photon, just as in e + e~ annihilation. In 
the cc and bb systems, assuming m a -C m c the branch- 
ing ratios are 

,«x M -(l^) 2 , (24) 
B (T(1S) -» 7 «) = (T(1S) -4 e*e~) 

We require that the contribution to J/^i — > 37 from 
J/ip — > 7a, a — )■ 77 is less than the experimental re- 
sult B (J/tp -> 3 7 ) = (1.2 ± 0.4) x 10" 5 07]. This trans- 
lates into a limit M > 6.2 GeV. Current limits on 
E(T(1S) -> 7a) at the 10" 6 to 10" 5 level assume a de- 
cays to leptons or hadrons [JS] and therefore are not con- 
straining on a pseudoscalar that decays promptly to 77. 
In addition, if a is long-lived enough to decay after pass- 
ing through the detector these decays are also subject 
to limits on J/tp, T(1S) — > 7 + inv as in the e + e~ — > 
7 + inv. case in Sec. |IVB| The current experimental 
limits are B {J/4> —> 7a — > 7 + inv.) < 4.6 x 10~ 6 for 
m a < 150 MeV gS] and B (T(1S) -> 7a -> 7 + inv.) < 
1.4 x 10 -5 for m a < 5 GeV [50] , To estimate the region 
of parameter space that these limits are sensitive to, we 
assume that the a decay length needs to be greater than 
1 m to go undetected in these searches. This excludes the 



regions 

/ 777 \ ^ 

from J/tp — > 7 + inv. and 

from T(1S) — > 7 + inv. 

We note that the regions of (m a ,M) ruled out by 
K + — » 7r + + inv and J/"0, T(1S) — > 7 + inv. are all con- 
tained within the region excluded by e + e~ — > 7 + inv. in 

Sec. EES 



D. Beam Dump Experiments 

Beam dump experiments, where weakly coupled parti- 
cles are searched for in the collision of proton or elec- 
tron beams with fixed targets, also provide limits on 
light pseudoscalars (often called axion-like particles in 
this context). To be successfully probed, the particles 
under consideration need to decay visibly and live long 
enough that they escape the target, but not so long that 
they do not decay before passing the detectors down- 
stream, typically tens of meters away from the target. 
This results in mass-dependent exclusion bands in the 
pseudoscalar's lifetime. 

Light pseudoscalars are produced in beam dump ex- 
periments either through the Primakoff process via their 
coupling to two photons or by being radiated via direct 
couplings to SM fermions as the beam constituents are 
stopped, in a manner analogous to bremsstrahlung. Each 
mode typically provides an 0(1) fraction of the total pro- 
duction cross section. The final states that are considered 
in searches for light pseudoscalars include 77, e + e~, and 

In the model we consider, the pseudoscalars are pre- 
dominantly produced by the Primakoff process since they 
do not directly couple to SM fermions and the dominant 
decay mode is 77. In placing limits on our parameter 
space, however, we simply apply the experiments' re- 
ported exclusions |51j to generic axion-like particles, ig- 
noring subtleties in the slightly different production cross 
sections and different branching fractions. This approxi- 
mation results in conservative exclusion regions since re- 
ducing the direct couplings to SM fermions should not 
make for more stringent limits. In the region of param- 
eter space of interest we find upper limits on the pseu- 
doscalar lifetime, corresponding to prompt decays inside 
the targets, or, equivalently, upper limits on the inverse 
coupling M in Eq. The values of these upper limits 
tend to roughly coincide with those coming from the re- 
quirement that the pseudoscalar also decays promptly at 
the LHC. 
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FIG. 5. Regions of pseudoscalar mass and inverse 077 cou- 
pling M excluded by the constraints in Sees. |IV"X|)IV E| The 
dashed curves are contours of 1 cm and 50 cm decay lengths 
for the pseudoscalar, assuming a boost 7 = m,h/2m a with 
m h = 125 GeV. 



E. Further Limits 

The coupling of a to SM fermions becomes negligi- 
ble in the non-relativistic limit. Therefore, the strict 
constraints light scalars and vectors are subject to from 
diffractive low energy neutron scattering on nuclei do not 
apply in this case. For the same reason, these particles 
are also not constrained by measurements of D-P tran- 
sitions in muonic Si and Mg. 

Additionally, light pseudoscalars can also be produced 
in the decays of excited nuclear states in nuclear reactors. 
However, this production mode is kinematically limited 
to m a < 10 MeV, and does not extend the exclusions we 
have considered. 



F. Allowed Regions 

In Fig. [5] we show regions of (M, m a ) excluded by the 
constraints outlined in Sees. |IV AlfTVE"] We also display 
contours of 1 cm and 50 cm decay lengths, assuming 7 = 
mh/2m a with m/, = 125 GeV. 



V. MODEL BUILDING 

Having investigated the phenomenology associated 
with a light pseudoscalar coupled to the Higgs, we spec- 
ulate on a few scenarios where the aFF interaction in 
Eq. |2]) could result. As we have seen, given a 125 GeV 
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Higgs, m a needs to be on the order of tens of MeV for the 
photons from its decay to be collimated enough to plau- 
sibly fake a single photon at the LHC. Using this fact in 
Eqs. ([5]) and |6]), we observe that M cannot be at the TeV 
scale for a to decay promptly as well to evade constraints 
from beam dump experiments. Obtaining a relatively 
low scale for a dimension-5 operator such as this involv- 
ing interactions beyond the SM poses a model-building 
challenge. 

A simple scenario where an interaction aFF is gen- 
erated is if a interacts with particles of mass m that 
have electric charge q in units of the electron's charge. 
One can imagine that a is either a composite made up of 
these charged particles or is fundamental and has a renor- 
malizable coupling to them. After integrating out these 
charged particles, the dimension-5 interaction between a 
and two photons is generated with a scale M ~ Air 2 m/q 2 . 
Given the requirement that a decay promptly seems to 
require a to couple to charged particles with a mass 
not much larger than tens of GeVjj Candidates for such 
charged particles are likely limited to SM fermions other 
than the t quark. Significant direct couplings to electrons 
or muons can be ruled out by anomalous magnetic mo- 
ment measurements. Couplings to heavy quarks c and b 
likely cause issues with well-studied quarkonium transi- 
tions. We are left with r and light quarks as SM media- 
tors of the aFF interaction. For a recent analysis of new 
light bosons decaying to 77 see Ref. [52] , 

1. Coupling to t 



limits g T < 2.5 for m a > 10 MeV. 

A similar limit of g T < 2 — 3 also comes from the OPAL 
Collaboration's result [M] 

B (Z ->■ t+t-j) < 7.3 x 1CT 4 , (31) 

where Z — > t + t~ a — >— > t + t~77 can mimic Z — > t + t~j 
events [S"2"] . 

Without a dedicated search for these pseudoscalars in 
r decays, it is difficult to be certain about constraints 
on this scenario coming from r decay. This is especially 
true if a has a significant invisible width, a scenario which 
we discuss in the Appendix. For example, the branching 
ratio of r to a and fi can be estimated to be |52j 

B (r ± -> a^vv) ~ 3.5 x KT 5 ^. (32) 

Some of these decays could appear as — > n^vwy de- 
cays, given the decay a — > 77. Even with the large 
Yukawa allowed by (g — 2) T and if all of the — > a^vv 
were registered as r — > fi Dvj decays, the contribution 
is still below the uncertainty on the measurement of this 
branching of (3.61 ± 0.38) x 10 -3 [55] which agrees well 
with expectations. Limits from r — > e ± Dvf are less 
stringent and hadronic r decays pose a bit more diffi- 
culty in prediction. 

In this scenario, the scale M is bounded by 

M = 47T ^- > 35 GeV. (33) 

9r 



We can also imagine that the effective aFF interaction 
is due to a coupling of a to t leptons, writing this coupling 



Car = -ig T O.TJ T. 



(28) 



This leads to an effective operator mediating a — ¥ 77 in 
Eq. (J2J with (assuming m a <C m T ) 



M 



A-n m T 

9r 



(29) 



The experimental information about the anomalous 
magnetic moment of r [5 3) 



-0.052 < ( g 2 ^ T < 0.013, 



(30) 



5 We note that a scenario where a — > 77 is mediated by a 
coupling of a to a new, light vector boson, V, that kineti- 
cally mixes with the photon (thereby lifting the requirement 
that a couple to electrically charged particles) is not viable 
because the amplitude for such a process is proportional to 

Q1Q2 (*?l — m v) 1 (^2 — m v) 1 wnere 91,2 are the photons' mo- 
menta, which vanishes for on-shell photons. If my vanishes, then 
a — > VV decays will strongly dominate the branching ratios. 

6 We write the interaction in terms of a pseudoscalar instead of an 
axial vector coupling here for simplicity. 



2. Mixing with ty° 

Coupling a to light quarks so that a mixes with 7r° 
could provide a solution. Given the couplings (ignoring 
s quarks for clarity) 



£ aq = -ig u auj u - ig d adj d, 



(34) 



using the partially conserved axial current, one can esti- 
mate the mixing angle, 



(9u-9d) \(qq)\ 
-FV (mlo - ml) 



(9u - 9d) 



(400 MeV) 2 



771 n — 777 1 



, (35) 



where (qq) is the value of the light quark condensate and 
F-n- = 92.2 MeV is the pion decay constant. Here, the 
scale M is given by its analogue in the w° case scaled by 
the mixing angle, 



M 



4tt 2 ^ 



sin ( 



(36) 



If there are no symmetry conditions on g u and gd, then 
we will in general expect additional isospin violation due 
to interactions involving a. The relative size of these 
effects is of the order g\ d j (^n 2 d /F 2 ^j . If we demand 
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that these not exceed ~ 1% then we require g Ut d < 10~ 3 . 
Furthermore, one must then contend with corrections to 
pion beta decay tt^ 1 



— > n e u which is shifted by 

T (tt* TT y, e ± v) = cos 2 6 r SM (tt ± 7T°e ± !/) . (37) 

The 0.6% agreement of theory and experiment for this 
rate [55] implies that 9 < 0.08. Such a bound is compat- 
ible with the limits above on the Yukawa couplings from 
isospin violation. 

The current measurement of the ir lifetime is at the 
3% level [57], and therefore does not pose any stricter 
constraints on the strength of this mixing. Furthermore, 
given a mixing angle of this size, the contribution to 7r° — > 



from a-7r° mixing is negligible. 



There are potentially strict limits in this scenario com- 



ing from the K + — > 7r + 77 limits described in Sec. IV C 



+ a is enhanced. The rate 



since the rate for K + — ) 
can be estimated as 

B (K+ -> TT+a) ~ sin 2 9 B (K + -> tt+tt ) (38) 

= sin 2 9 x 0.21. 

The limits on this branching at the 10 -7 level imply 
that 6 < 10~ 3 , or equivalently, M > 1 TeV. This does 
not allow for prompt decays at the LHC if a is produced 
in the decay of a 125 GeV Higgs and is problematic with 
respect to beam dump constraints. However, the anal- 
ysis in Ref. [IS] only selected events corresponding to 
m a < 100 MeV because of the very large background 
from K + 7T + 7T . Thus, to avoid the strict limit on 
M that results from this analysis, a must have a mass 
greater than 100 MeV if its coupling to photons is medi- 
ated by a mixing with ir°. 

The limit from pion beta decay implies that M > 
47 GeV in this case. 



3. Other Scenarios 



rate depends sensitively on the photon identification cri- 
teria used in the analysis. Furthermore, the presence of 
more photons implies that the fraction of events with a 
photon conversion is higher than in a pure 77 sample. Al- 
though CMS has already presented the best-fit rates in 
conversion and unconverted events separately, the statis- 
tics remain low, and we leave such analysis for the future. 
We also note that in the case where the a is long-lived, 
7cr ~ 50 cm, it is imaginable that the decay length may 
be determined from the distribution of conversion radii, 
even though since the states are highly boosted they do 
not display displaced vertices in the traditional sense. 

The light pseudoscalar hypothesis can also be tested 
by looking for a 37 signal with invariant mass matching 
that of the Higgs boson. If the Higgs is produced with 
a sizable boost, some decays will yield one slow a and 
one fast a in the lab frame, and the former will produce 
widely separated photons. Alternatively, if m a is much 
larger than the values studied here, the pseudoscalar can 
still contribute some events to h — > 27 if B(h — > aa) is 
large enough, but then the particle may be more easily 
seen in searches for h — > 47 |16) . 

A different way to search for a light pseudoscalar cou- 
pling to photons is to produce the a directly in Primakoff- 
type experiments, 7+Nuc. — > a+Nuc, which would allow 
the mass of a to be directly reconstructed. A proposal 
for an upgrade to the PrimEx Experiment at Jefferson 
Lab to measure T (r) — > 77) envisions collisions of pho- 
tons with an energy E 7 ~ 11 GeV on a liquid hydrogen 
target with a luminosity of C ~ 10~ 2 nb _1 s _1 with a run 
of 45 days [55] • Using the interaction in Eq. ^ , the cross 
section for the Primakoff production of a in the collision 
of a photon with a proton can be found to be 



da _ 2^ 2 a 3 |F cm (Q 2 
dQ 2 ~ M 2 Q 2 



(39) 



In addition to the two scenarios described above we 
briefly mention a few other potential models that could 
lead to the interactions in Eq. ([2| with a scale M on 
the order of tens to a few hundred GeV. As in [T3], the 
pseudoscalar could be coupled to heavy vector-like mat- 
ter that mediates its interaction with photons. However, 
ensuring that the scale M is small enough in this situa- 
tion so that a decays promptly is difficult. This could be 
solved by giving the heavy matter a large electric charge 
or gauging it under a non-Abelian group with a large 
number of colors. We also speculate that sterile neutri- 
nos with very large transition magnetic moments coupled 
to a could also offer a solution since their masses need not 
be at the weak scale. 

VI. OUTLOOK 

As discussed in Sec. |IIIB] the scenario analyzed in this 
paper can be tested in h — > 77 events, since the expected 



where Q is the momentum transferred to the proton and 
F em is its electromagnetic form factor. Using a simple 
dipole form factor and a proton charge radius of 0.87 fm 
(the results do not depend sensitively on the form factor 
and charge radius — using 0.84 fm as measured in muonic 
hydrogen does not change the estimate) and a detection 
efficiency of 60%, this cross section would lead to 

a AO GeV\ 2 , x 

^10^— j (40) 

pseudoscalars collected with m a = 40 MeV. For com- 
parison, there would be about 10 ?7's produced via the 
Primakoff process and subsequently decaying to two pho- 
tons under the same conditions. 

It appears likely, at least statistically, that pseu- 
doscalars with parameters in the range that is interesting 
in the context of the h — ¥ 77 signal can be probed at fu- 
ture Primakoff experiments. 



11 



VII. CONCLUSIONS 



Appendix: Invisible Decays 



In this paper we have investigated a simple model that 
can give rise to an apparent excess relative to the SM in 
the h — > 77 channel. In our study the two photon ex- 
cess comes from Higgs decays to two light, boosted pseu- 
doscalars, h — » aa, and each pseudoscalar decays into 
two photons. For very light pseudoscalars, each pair of 
photons is highly collimated and a non-negligible fraction 
of the pairs can appear as single photons, even at high- 
resolution detectors like those at the LHC. This scenario 
serves as an example where photon jets [T5] are produced 
and helps to motivate the future experimental study of 
these objects in more detail. 

We have estimated the fraction of events that ap- 
pear as 27, taking into account the fine granularity of 
the ATLAS detector and the tight photon selection cri- 
teria. Additionally, we have investigated subtleties that 
may arise when one photon or more converts. We have 
assumed that the different analyses used at CMS and 
ATLAS do not give rise to large deviations between the 
h — > 77 measurements, but we note that the amount of 
47 contamination can be analysis-dependent and thus the 
differing analyses may be used as a test of the hypothesis. 

In addition to a diphoton excess, the model also pre- 
dicts a deficit relative to the SM for Higgs decays into 
other final states. Although the uncertainties are large, 
the data that have been collected so far at the LHC and 
Tevatron are consistent with both of these predictions, 
and with more data the model will be easily tested. The 
relative size of a potential excess in the 77 channel and 
a decrement in the remaining channels will be a direct 
probe of the mass of the pseudoscalar, with a lighter 
pseudoscalar better able to accommodate a larger signal 
in the non-77 channels for a fixed 77 rate. Lower bounds 
around m a ~ 10 MeV from 7 + inv. searches and beam 
dump experiments offer a complimentary sensitivity to 
the model. 

Finally, the mass of the light pseudoscalar can be well- 
measured at Primakoff experiments. As we discussed in 
Sec. |VI| future Primakoff experiments appear well-poised 
to test models involving new light bosons coupled to pho- 
tons like the one we have considered, underscoring the 
complementarity of experiments at the intensity frontier 
with high-energy studies of the Higgs boson. 



In this Appendix, we describe a simple extension of 
the model in Sec. |ll] to allow for a to decay invisibly. 
The Higgs phenomenology and direct constraints change 
somewhat from the case where a decays purely to 77. 
As before, we present fits to the Higgs data and show 
allowed regions of parameter space. 



1 . Framework 

In addition to the two photon decay channel, the 
pseudoscalar a could have an appreciable invisible decay 
width. As a concrete example, we can add an interac- 
tion with a stable (on collider scales) SM singlet Dirac 
fermion x, 



r 



M 



7 X7 M 7 5 X, 



(A.l) 



where M' is another scale describing the strength of this 
interaction. The rate for a — > \X 1S then 



mZm n 
r (a -> XX) = 7^P X 



When this channel is open, the branching ratio for a 
77 becomes 

0-1 -1 



(A.2) 



B (a — > 77) = 



1 



p x ( 1 m x M 
■kol m„ M' 



(A.3) 



where pi = 1 — 4m^/m„. The decay length (for a boost 
7 = m,h/2m a ) is shortened by a factor of this branching 
ratio, 

7 cr ~= 1.15 mm ( T77 ^ T7 ) ' f A ™°^ T ) ^ (A.4) 



10 GeV/ V 40 MeV, 

( m h \ ( B (0 -> 77) \ 
U25 GeV/ V 0- 1 / 



x 



.125 GeV, 

and the scale M can then be expressed as 
, 1 / 2 / m„ 



M = 29GeV (J^LY" ( 
VI cm/ V 



-V 



40 MeV/ 
m h \-Va / g(a-> 
125 GeV/ V 0.1 



(A.5) 

77)r 1/2 . 



For a fixed decay length, a sizable invisible branching 
fraction for a allows the scale M to be somewhat larger 
than it would be in the case where a only decays to pho- 
tons. 
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2. Higgs Signal 

The primary changes to Higgs phenomenolo; 
B (h — > 77) eff and i? 77 , now given by 

B i h -> 77) e fi = B(h^ 77) 



igy are in 



- eB(h—^ aa) B (a —> j-f) 2 



(A.6) 
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FIG. 6. Regions of light pseudoscalar parameter space that 
are favored (blue/light) and disfavored (red/dark) by the cur- 
rent best-fit signal strengths in the h — > 77, ZZ, WW, bb, tt 
channels. Contours are overlaid for the net diphoton (solid 
green) and ZZ, WW, bb, tt rates (dashed yellow) expected at 
the LHC relative to the SM rates. B (a — > 77) is fixed to 0.1, 
indicating a 90% branching of a — > invisible. 
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FIG. 7. Same as Fig. [5] in the case that B (a — > 77) = 0.1. 
The dashed curves are contours of 10 and 50 cm decay lengths 
for 7 = mh/2m a and mi, = 125 GeV. To the left of the 



vertical dotted line, R 11 > 1 [Eqs. |7| and (111] while to 



-R 77 < 1. In this case, only the beam dump and 
7 + inv. limits, which are the most constraining, are 



shown. 



and 



Direct Constraints on a 



1 + B (/!->■ aa) 



eB (a — y 77) 
Bsm (h -> 77) 



(A.7) 



Fixing B (a — > 77) = 0.1, we regenerate the exclusion 
contours and contours of fixed R~ n ,xx in Fig. [6] The 
best-fit x 2 /d.o.f. = 6.5/10, and the best-fit point is close 
to the intersection of the i? 77 = 2 and Rxx — 0.5 con- 
tours. Like the B (a — > 77) = 1 case, the \ 2 is shallow 
and the best-fit point is not yet a meaningful quantity, 
so we do not show it in Fig. [6j 

In addition to an increased apparent branching to 77, 
there are contributions to the invisible width of the Higgs, 
which can be related to the increased diphoton branching 
ratio, 

B (h -> invisible) = B (h -> aa) B (a -> xxf (A.8) 
= B{h^aa) [1 - B (a -> 77 )] 2 . 

There is also a "monophoton" branching, 

B (h -> 7 + invisible) = B (h -> aa) B {a -> 77) (A.9) 

x B (a -> XX) x 2e' 
= 2e' B(h^ aa) B (a -> 77) [1 - B (a -> 77)] , (A.10) 

where e' is the efficiency for the one photon jet in the 
event to be reconstructed as a single photon. Future 
LHC monophoton searches with low MET cuts may be 
able to probe this channel. 



We now analyze changes to the direct limits given in 
Sees. |IV A] |IV E| that arise if a has an appreciable invisible 
branching fraction. 

Signals involving missing energy now apply regardless 
of a's decay length. The e + e~ — > 7 + inv. limit then 
becomes the strongest lower limit on M. 

The limits on M from beam dumps must be rescaled 
by a factor of B (a —> 77) since they are limits on 
the a lifetime. We conservatively ignore any shrinking of 
the exclusion region due to the smaller visible branching 
ratio. 

An invisibly decaying a could contribute to star cooling 
by providing an additional channel for energy to leave a 
star. However, such limits are not important in the mass 
range that we consider here. 

If a mixes with the 7r° and can decay invisibly, there are 
stringent limits on the invisible decay rate of the 7r°, with 
the collider upper limit measur cd to be 2.7 x 10~ 7 [59] 
and a limit from Big Bang Nucleosynthesis that is several 
orders of magnitude stronger |60j . These are quite con- 
straining and make the scenario where the aFF cou pling 
is generated by a mixing with 7r° described in Sec. V 2 
unlikely if a has an appreciable invisible branching frac- 
tion. 

A large invisible branching fraction for a renders the 
model where the coupling to photons is mediated by a 
coupling to t more plausible. This is because r decays 
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necessarily involve at least one neutrino in the final state We reproduce the low-energy exclusions of Fig. [5] in 
so observing missing energy in a r decay does not signal Fig.[7]for B (a — > 77) = 0.1, showing only the most strin- 
new physics or a rare SM decay. gent limits from beam dumps and e + e~ — > 7 + inv. 
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